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Abstract 

In this paper, we generalized the surface integral equation method for the evaluation of Casimir 
force in arbitrary three-dimensional geometries. Similar to the two-dimensional case, the evalua- 
tion of the mean Maxwell stress tensor is cast into solving a series of three-dimensional scattering 
problems. The formulation and solution of the three-dimensional scattering problem is well-studied 
in classical computational electromagnetics. This paper demonstrates that this quantum electrody- 
namic phenomena can be studied using the knowledge and techniques of classical electrodynamics. 
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The Casimir force has first been predicted by Casimir in 1948 jlj. It exists among 
charge-neutral bodies due to the quantum fiuctuation of the elecromagnetic fields in vacuum. 
Although it is tiny when the objects have large separations, it becomes the dominant force 
between charge-neutral objects when the separation is below micron meters. 

Micro-electromechanical systems (MEMS) are micron-sized devices in which mechanical 
elements and moving parts, such as tiny sensors and actuators are carved into a silicon 
substrate. As further miniaturization takes place, the devices may reduce to nano-scale and 
become nanoelectromechanical systems (NEMS). They have a wide class of applications. For 
example, the release of the airbag in cars is controlled by a MEMS-based accelerometer. One 
of the principal causes of malfunctioning in MEMS is stiction, i.e., the collapse of movable 
elements into nearby surfaces, resulting in their permanent adhesion. Casimir eflFect is often 
an important underlying mechanism causing this phenomenon [2^. On the other hand, this 
eflFect could be put into good use. Capasso and his group showed that the force can be used 
to control the actuation of a micromachined torsional device jsl . This group also showed that 
the Casimir attraction can be used to make a nonlinear oscillator. The force infiuences the 
dynamical properties of a micromachined device, changing its resonance frequency, hysteretic 
behavior, and bistability in its frequency response to an AC excitation. They proposed that 
this device could serve as a position sensor working in nanometer scale jj, [5]. 

Since Casimir effect is rather significant in MEMS and NEMS devices, it is important to 
be able to predict it accurately before the devices going into manufacturing. The geometry 
of the MEMS and NEMS devices may be complicated; thus approximation methods, such 
as the proximity force approximation (PFA) and optical approach in terms of virtual 
photons moving along ray optical paths [7], are no longer sufficient. Instead, an exact and 
general numerical method that can handle arbitrary geometry is desired. With the numerical 
method, we can get the desired accuracy by giving enough computational resource. Recently, 
major progress has been made in the development of numerical methods for Casimir force. 
In general, there are two popular approaches: the path integral approach (also known as 
EGJK method) and the Maxwell stress-tensor approach. 

Schwinger is the first one to attribute the Casimir interaction to fiuctuating current 
and charge densities inside the objects. With path integral method, the Casimir energy 
is expressed as the logarithm of a functional integral over all allowed configurations of the 
fiuctuating currents, weighted by the appropriate action [8]. It has been generalized in [9] 
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(EGJK) for compact objects of arbitrary shape and separation and applied to predict the 
force between a cyhnder and a plate [l^, and between two spheres j^. EGJK's agorithm is 
efficient to predict the Casimir energy between compact 3D objects of spheroidal or nearly- 
spheroidal shape. The algorithm leads to analytically tractable series solution. But it is 
not of practical use for general geometries due to the poor convergence rate for objects 
with corners. This limitation occurs since they used spherical basis functions to expand the 
current distribution. It could be overcome by expanding the currents with basis functions 
defined on a pair of adjacent triangles, which is known as the RWG basis {ll,[l2]. The RWG 
basis is vastly popular in the surface integral method in computational electromagnetics. 
After some manipulation, the evaluation of Casimir energy depends only on solving for all 
the eigenvalues of the impedance matrix obtained from the integral equation by the method 
of moments. 

The theoretical foundation for the Maxwell stress tensor approach was built by Lifshitz 
and his coworkers 1^: the net Casimir force on a body can be expressed as an integral 
over any closed surface around the body of the mean electromagnetic stress tensor, inte- 
grated over all frequencies. And the fiuctuation-dissipation theorem states that the mean 
fiuctuating field, which forms the Maxwell stress tensor, is related to the imaginary part 
of the dyadic Green's function of the same geometry. The dyadic Green's function can be 
evaluated by existing numerical methods in computational electromagnetics (electrodynam- 
ics) {14]. Rodriguez has demonstrated the applicability of this method by using a simple 
finite-difference frequency-domain (FDFD) method to calculate the dyadic Green's func- 



tion in z-invariant structures 



15 



16| . This approach relates the quantum electrodynamic 



phenomena with classical electromagnetics which motivates our piece of work: to evaluate 
the Casimir force among arbitrary objects using integral equation method. With surface 
integral equation method, the number of unknowns can be greatly reduced. Moreover, fast 
algorithms [17] are available to further improve the efficiency of the method. The compu- 
tational complexity of both the stress tensor approach and the path integral approach can 
be reduced to 0{N log N) times the number of iterations at best, being the number of 
surface unknowns. Thus their efficiency is comparable. The main difference between them is 
the physical quantity being calculated directly: path integral method calculates the energy 
while stress tensor approach calculates the force distribution (pressure). 

In a paper published earlier [20], we have demonstrated that surface integral equation 
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method could be applied to evaluate the dyadic Green's function and its derivatives in 
two-dimensional (2D) geometries. By using this method, the Casimir force between 2D 
objects could be simulated using much fewer unknowns and can be easily accelerated by 
fast-algorithms readily developed in computational electromagnetics. In this paper, we 
would extend the surface integral equation method for arbitrary three-dimensional (3D) 
geometries. The idea is similar to the last paper |2o|: the Casimir force is expressed in 
terms of the dyadic Green's function and its derivatives and they could be calculated from 
the integral equation by using different types of dipole point sources. 

The starting point of the stress tensor approach is that the net Casimir force actin g on 



the surface S of an object is given by an surface integral of the mean stress tensor on it 

F = ^(T(rO) -rfs' (1) 



where T is the Maxwell stress tensor defined as in [18]. For an arbitrary 3D object, we 
define u and v as the tangent vector at the surface, and h as the normal direction of the 
surface. They satisfy the relationship it x ij = h. If the object is a perfect conductor, the 
tangential electric field Ey and the normal magnetic field Bn vanishes at the surface of 
the object. The stress tensor could be simplified as follows: 

Tnnir) = leoElir) - -^5^(r) - -^5,^(r) (2) 

Tun{r) = r,,(r) = (3) 

Thus for perfect conductors, the Casimir force becomes a pressure on the surface. 

The average of the fiuctuating electric and magnetic fields in the ground state is obtained 
from the fiuctuation-dissipation theorem fiol : 

(0|£;,(r,^)£;,(r^^)|0) = -Im / c^2Q,(r,r^c^)rfu; (4) 

(0|5,(r,^)5,(r^^)|0) = -Im / (Vx),,(Vx),^Q^(r, r^ c^)dc^ (5) 

^ Jo 

where the dyadic Green's function is the same as that defined in classical electromagnetics 

n 

to relate the current to the field in an arbitrary geometry [20]. 

In 3D structures, the field components and the dyadic Green's function satisfy the fol- 
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lowing relationship: 



E^{r) = lu [ dv'G^p{k,r,r')Jp{r') 
Js 

Huir) = - [ dr' 
/^o Js 



d d 



(6) 
(7) 



Substituting (|6]) and ([7]) into (|4]) and (Ej) , the average of each term in the stress tensor could 
be represented by 



-eo(^.(r)E„(r')) 



1 
1 

2/io 



hco ^ r ki 

Im / —En{ko,v,Si)dkQ 

^Im / dko[-Hu{ko,r,S2) + Hy_{ko,r,S3)] 
[B^{r)B^{r')) = ^Im dko[-H^{ko,r, S4) + Hy_{ko,r, s^)] 



[B^{v)Bu{r')) 



(8) 
(9) 
(10) 



Here, E(/co, r, Si) and H(/co, r, Si) are the electric/magnetic field at location r in the geometry 
under the excitation of source Si. Their corresponding relationships are summarized in Table 



TABLE I: A table of source and field operators for each term in the stress tensor. 





Field Component 


Source (r^ = r) 


(1) 


En{v,Si) 


hS{r' - Ti) 


(2) 




V irSiv' - r,) 


(3) 


Hu{r, S3) 


n iyS{v' - r,) 


(4) 


i/^(r,S4) 




(5) 







If the object is made of perfect conducting material, there is only induced electric current 
J fiowing on the surface. The scattered electric field and magnetic field in the free space 
region could be obtained from the induced electric current by C and JC operators respectively 
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E(r) = E*(r) + /:j(r) 
H(r) = (r) - 0(r) 



(11) 
(12) 



where E and H are the total field in the space and and are the incident field generated 
by an external source. The operators are defined as 



£X(r) 



1 + 



w 



g{r,r')-:^{r')dr' 



/CX(r) = J^X(r') X Vg{r,r')dr' 



(13) 
(14) 



where ^'(r, r') = /47r|r — r'|. The total tangential electric field vanishes on the surface 

of the PEC object. Enforcing Eqn. ( |T3l ) on the surface we obtain the electric field integral 
equation(EFIE): 



h X E'(r) = h X C3{r) 



(15) 



By expanding the unknown current with RWG basis jl2], J(r') = X^n ^"*'"(''') ^^*1 
ing Galerkin's method, the integral equation ( fT^ could be cast into a matrix equation 
[Zmn\[an\ = [9m], whcre 



drJ^(r)-EXr) 



Zry 



' Tm 





= vv" 


dr dr'J^ ■ 


^0 _ 


J Tm ^ Tn 



^(r,rO-Jn(r') 



(16) 
(17) 



If the sources listed in Table [J are used to generate the incident field, then the evaluation 
of involves integrals containing some super-hyper singularities. The treatment of these 
singularities has been discussed in detail in j22|]. 

Now we would present some numerical results obtained from this method. They are 
compared with either analytical results or other's calculation. First we would start with the 
evaluation of the Casimir force between two spheres, since they are the simplest 3D objects 
and the analytical expression of the Casimir energy is available in [q]. We assume that the 
radius of both spheres is R and their minimum distance is Z, the product of the force and 
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the square of the radius (F* R'^/hco) is a dimensionless and scale-invariant quantity. So the 
result is given by their product (F^R^/hco) versus the normalized separation (Z/R). In the 
numerical process, first the Casimir force per unit area (pressure) is calculated at different 
points on the surface, and then the total force is obtained by integrating the pressure over 
the entire surface. The analytical results are interpolated from the data given in ref. 
Good agreement has been observed. 
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FIG. 1: Casimir force between two spheres of radius R and a separation of Z. 

The geometry of the second example contains two identical capsules. A capsule is a cylin- 
der of radius R with hemispherical end caps. Its total length is denoted by L. Experiment 
had been carried out to measure the interaction between two crossed cylindrical capsules (23 1 
so a numerical study of the geometry would be interesting to experimentalist. It has also 



been simulated by using EGJK method in ref. [UJ] and reference data is available to make 
comparison. The force between two capsules is evaluated for two cases: they are parallel or 
they are perpendicular. ^ ^ 

Similar to the 2D geometries [20] case, the integral equation method could reduced the 
computation cost significantly since it involves only surface unknowns. Moreover, both the 
coordinate space integration ([1]) and the spectral space integration (|8j) to (fTOl) are smooth 
and independent of the number of unknowns. 

As we know, both our surface integral method following the stress tensor approach and 
Reid et al numerical method following the path integral approach [11] are eflFective numerical 
methods for calculating the Casimir force among arbitrary 3D objects. Though they start 
from completely diflFerent description of the Casimir energy and force, they are connected 
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FIG. 2: Casimir for two pairs of perfectly conducting objects: identical capsules with parallel axes, 
and identical capsules with perpendicular axes. Here L — QR. SIE refers to the surface integral 
equation method derived in this paper. 

with the surface integral equation method in classical electromagnetics at certain stage. 
Reid et al method reduces the problem to a search of the all eigenvalues of the impedance 
matrix. Our method reduces the problem to solving a set of matrix equation, with the same 
impedance matrix as in Reid et al method. Reid et al method calculates energy directly 
while our method calculates the force distribution directly. The complexity of both method 
could be reduced to 0(A^log A^) at best by using fast algorithms, where is the number of 
surface unknowns. The drawback of the path integral method is that we need to differentiate 
the total energy with respect to displacement in a certain direction to get the force. This 
operation changes the problem into a generalized eigenvalue problem which increases the 
cost and complexity. Another drawback is that the force distribution on the surface is not 
available from the method. On the contrary, for stress tensor approach, force distribution is 
required before we obtain the total force. The tradeoff is that an additional surface integral 
needs to be performed compared to the path integral method. 
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